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"c/2 Abstract 

jS^ During therapeutic treatment with heavier ions like carbon, the beam under- 

I ^-h goes nuclear fragmentation and secondary light charged particles, in particu- 

lar protons and alpha particles, are produced. To estimate the dose deposited 
^ into the tumors and the surrounding healthy tissues, the accuracy must be 

cn higher than (^±3% and ± 1 mm). Therefore, measurements are performed to 

determine the double differential cross section for different reactions. In this 
paper, the analysis of data from ^"^ C +^'^ C reactions at 95 Me V/u are pre- 

•^ sented. The emitted particles are detected with AEthin-^Ethick-E telescopes 

made of a stack of two silicon detectors and a Csl crystal. Two different 
methods are used to identify the particles. One is based on graphical cuts 
onto the AE-E maps, the second is based on the so-called KaliVeda method 
using afunctional description of AE versus E. The results of the two meth- 
ods will be presented in this paper as well as the comparison between both. 
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1. Introduction 

Carbon ions witli energies ranging from 80 to 400 MeV/u are used in 
hadrontherapy to treat cancerous tumors. Their advantage compared to 
photons hes in the great accuracy of ion beams to target the tumor while 
sparing the surrounding healthy tissues (due to the high dose deposition in 
the Bragg peak and the small angular scattering of ions). The physical dose 
deposition is affected by the fragmentation of the ions along their penetration 
path in the patient tissues [T]. This process leads to a consumption of the 
primary beam (up to 70% for 400 MeV/u ^^C in water) and to a delocalisation 
of the dose in the healthy tissues |2j . 

Simulation codes are used to compute the transport of ions in matter, 
but the constraints on nuclear models and fragmentation cross sections in 
the energy range used in hadrontherapy (80 to 400 MeV/u) are not yet suf- 
ficient to reproduce the fragmentation processes with the required accuracy 
for clinical treatments. To improve the knowledge on the ^^C fragmentation 
process, experiments have been performed in Japan and in Europe for more 
than 15 years. Measurements of light charged fragment production in water 
and PMMA (C5H8O2) have been started by the Japanese treatment centers 
(Chiba and Hyogo) and by the Gesellschaft fiir Schwerlonenforschung (GSI) 
biophysics department in the energy range 200-400 MeV/u [3] [1] [5] [H] [Z] [H] ■ 

To extend these data to the lowest energies, a first integral experiment 
on thick water equivalent targets has been performed by our collaboration 
in May 2008 at the Grand Accelerateur National d'lons Lourds (GANIL, 
France). Energy and angular distributions of the fragments produced by 
nuclear reactions of 95 MeV/u ^^C with thick PMMA targets have been ob- 
tained p]. Comparisons between experimental data and simulations, using 
different physics processes of GEANT4 |10] for modeling the passage of parti- 
cles through matter, have shown discrepancies. Nevertheless, it is difficult to 
constrain the nuclear reaction models by a direct comparison to thick targets 
experiments. 

To improve the models and reach the accuracy required for a reference 
simulation code for hadrontherapy, a second experiment has been performed 
on thin targets on May 2011 at GANIL to study C-C, C-H, C-0, C-Al and 
C-""*Ti reactions at 95 MeV/u. The particle identification has been per- 
formed by using telescopes (stack of detectors) and the AE-E technique. 
The experimental set-up will be presented in the next section. 

In such experiments, the usual procedure to extract the particle natures 



and the particle energies from raw data consists in two separate steps: the 
particle identificaton and the energy calibration of the detectors. The particle 
identification is done by using graphical cuts on so called AE-E maps. The 
energy calibration is done by using well defined reference energies (elastic 
energies, punch through energies). We propose in this paper an alternative 
procedure which performs the particle identification and the detector calibra- 
tions in one single step. The two different analysis methods will be described, 
the first one called "graphical method" uses graphical cuts, the second one 
called "KaliVeda method" uses the KaliVeda framework developed by the 
INDRA collaboration [TT] [12] . The main advantage of the KaliVeda method 
is that once the identification grids are created, they allow to do both the 
detector calibration and the particle identification at the same time thanks 
to automated procedures which saves time in data analysis. A comparison 
of the results obtained with the two methods for the C-C reaction will be 
presented in the last section. 

2. Experimental setup 

2.1. Charged particles detection 

The experiment has been performed using the ECLAN reaction chamber 
at the GANIL G22 beam line. The energy of the impinging carbon ions was 
95 MeV/u on different thin targets (C, CH2, Al, AI2O3, ''"^Ti) with area 
densities of about ~100 mg/cm^. In this section we will report only about 
carbon target data. The setup consists of four AEthin-AEthick-E telescopes 
mounted two by two on two rotating stages that allows rotation inside the 
chamber from 7° to 43° at a fixed distance of 22 cm behind the target. Thus 
the solid angle covered by each telescope was 6.5 msr. A fifth telescope was 
at fixed angle at 4°, downstream, covering a solid angle of 0.4 msr as shown 
in Fig. [Tj 

The five AEthm-AEthick-E telescopes consist in two silicon detectors fol- 
lowed by a 10 cm thick Csl scintillator. The thicknesses of the silicon detec- 
tors were 150 yum and 1000 /xm respectively. 




(a) (b) 

Figure 1: Representation of the experimental setup. 



2.2. Beam Monitor 

The beam monitoring is based on the measurement of fluorescence X- 
rays emitted by a thin Ag foil (7 /xm thick) set in the beam after the reaction 
chamber. X-rays are detected by means of two Si (Li) detectors located at 
90° with respect to the beam direction. These two detectors are calibrated 
at different beam intensities ranging from 10^ to 10^ ions/s by means of a 
plastic scintillator intercepting the beam. The beam intensity is determined 
with the two detectors (Si(Li)). The detected fluorescence X-rays, emitted 
by the thin Ag foil, are represented by two peaks in the energy spectra 
of the Si(Li) (Fig. 2(a)). These peaks have been integrated while taking 
into account the background. Different background subtraction methods 
have been tested: no background, constant background, linear background, 
polynomial background. They are all compatible within 2%. The energy 
spectrum of the plastic scintillator presents several peaks corresponding to 
the number of carbon ions detected in coincidence (Fig. 2(b)). These peaks 
have been fitted with gaussian functions and integrated hereafter. Fig [3] 
represents the beam intensity in the plastic detector relative to the detected 
X-rays intensity in one Si (Li) for several calibration runs. The calibration 
parameters are obtained by a linear fit of the calibration points. This monitor 
gives a global accuracy of 5% on the whole intensity range of the experiment 
(from 10^ to 10^ ions/s). 
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Figure 2: Left-hand side: Energy spectrum of the X-rays emitted by the Ag foil obtained 
with a Si(Li) detector. Right-hand side: Energy spectrum obtained with the plastic scin- 
tillator for the beam monitor calibration. 
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Figure 3: Intensity calibration of the SiLi detector thanks to the plastic scintillator inter- 
cepting the beam at low energy (Ipiast — pOxIgi(Lj) -f pi). 



3. Data Analysis using graphical cuts 

3.1. Particle identification 

The AE-E technique is used to perform the identification of the charged 
particles. This method is based on the fact that the energy loss of charged 
particles in matter depends on their charge and mass. Figs. |4] and [5] show 
examples of identification maps for the 250 /im thin Carbon target. Each 



hyperbole corresponds to a particle characterized by its charge and mass 
numbers (Z,A). Figure |4] represents the energy deposited in the thin silicon 
detector relative to the energy deposited in the thick one. The hyperbolas 
back bending are due to the particles which cross the second stage of the tele- 
scope. Figure |5] shows the energy deposited in the thick Si detector relative 
to the energy deposited in the Csl scintillator. This map is complementary to 
the previous one, thanks to the identification of the particles passing through 
the first two stages. The identification is performed by a set of contours for 
each telescope in order to separate each particle in charge and mass, from Z 
= 1 to Z = 6. These contours are represented by black lines on the maps 
of the figures |4] and |5| Since the data acquisition is triggered when the sec- 
ond stage is hit, the identification threshold at low energy is determined by 
the thickness of the first stage. The minimum energy requisite to cross the 
first stage (150 fim.) and thus to hit the second one depends on the nuclear 
charge of the particle, these threshold values are given in Table [T] The Csl 
crystal thickness prevents most of the charged particles from passing through 
the scintillator. Indeed, the lightest particles namely protons, deuterons and 
tritons, need respectively energies greater than 198, 130 and 103 MeV/u to 
cross the Csl crystal. 



Nuclear charge Threshold (MeV) 

i 4 

2 14 

3 30 

4 46 

5 61 
6 81 

Table 1: Minimum energy for detection as a function of the charge of the particles. 
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Figure 4: Identification map between the first stage and the second stage of a telescope. 
The black lines represent the different graphical contours used for the particles identifica- 
tion. 
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Figure 5: Identification map between the second stage and the third stage of a telescope. 



3.2. Energy calibration 

The energy calibration of the two sihcon detectors is based on the maximal 
deposited energy in the second stage for a given particle without crossing it. 
Thanks to the known thicknesses of the Si detectors, this energy threshold 



is computed, for each particle, by using SRIM [13] simulations and is linked 
to the corresponding channel number of the thick Si detector. The energy 
loss in the thin Si detector is deduced from the energy deposited in the thick 
Si detector. This allows to establish the correspondence between the energy 
deposited in the thin Si detector and the channel number in the corresponding 
charge to digital converter (QDC) module. Fig. M shows an example of the 
calibration points obtained for one thick Si owing to this method. As the 
response of a silicon detector does not depend on the nature of the particle, 
the correlation between the energy loss and the channel number is the same 
for all particles. The calibration parameters are obtained by a linear fit of 
the points as shown on Fig. |6} 

Contrary to silicon detectors, the produced light in the scintillator de- 
pends on the nature of the particle and on its energy. It is necessary to 
use a calibration function with a distinct set of parameters for each isotope. 
The chosen function, based on the one proposed by Stracener et al [H], is 
described by the following formula: 

EmeV = a + /5Echannel + 'yln{l + (5Echannel) + ^Echannel (l) 

This equation takes into account the quenching of the scintillator at small 
energy as well as the linear response at higher energy. The term, rjE^^^^^^^, 
has been added to take into account a possible saturation at very high energy. 
The energy loss in the Csl detector is deduced from the energy deposited in 
the thick Si detector using SRIM simulations. This method propagates the 
errors made on the energy calibration of the thick Si detector to the energy 
calibration of the Csl detector. The lower the deposited energy in the thick 
Silicon detector is, the larger the error on the residual energy released in 
the Csl is. With our experimental setup, this error has been estimated, 
from almost 30% for 95 MeV protons to 4% for 95 MeV/u carbon ions. An 
example of the calibration functions obtained is shown on Fig. [7j The five 
parameters (a, (3, 7, 6 and rj) are determined for each isotope and for each 
telescope. 
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Figure 6: Energy calibration for the thick sihcon detector of a telescope (linear fit). 
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Figure 7: Calibration of the Csl scintillator of a telescope thanks to the calibration function 
given by Eq. [T] 



4. Data Analysis using KaliVeda 

An new analysis method based on the KaliVeda toolkit |12] has been used 
for the analysis of the experiment. KaliVeda is an object-oriented software, 
based on the ROOT [TJ] framework for data processing, with many useful 
classes for data analysis and energy loss calculation. It has been developed 



9 



to analyse data obtained by the AE-E technique. It was firstly developed for 
the analysis of the INDRA multi-detector. 

4.1. KaliVeda tools 

In this part, the KaliVeda tools which are used for the analysis will be 
presented. 

4-1.1. The Energy Loss functional 

KaliVeda is using a functional describing accurately the energy loss AE 
in a detector as a function of the residual energy E deposited in a second 
detector in which the particle has stopped. This functional is the one pro- 
posed by L.Tassan-Got in which the charges Z and the masses A of particles 
are parameters [16]. For detectors delivering a linear response (as silicon 
detectors), it reads: 

AE = [(^E)'^+^+i + (AZ"A^)'^+'^+i + ^Z2A'^(^E)'^] ^^+^ - ^E (2) 

where g, /i, v, A, a, /3, ^ are fit parameters. 

For detector delivering a non linear response versus the deposited energy 
this functional needs to be corrected. For scintillators, it allows to take into 
account the quenching effect of scintillation leading to a non linearity of the 
emitted light versus the deposited energy. The energy E in equation [2] must 
now be expressed as a function of the light h emitted by the scintillator. The 
light response of CsI(Tl) crystals is deduced from the Birks formula [T7] and 
allows to express the energy released in the Csl as a function of the emitted 
light : 



E= \lh^ + 2ph 
where p = rfl^K is a new fit parameter. 



1 + /n ( 1 + - 



(3) 



4.1.2. The KaliVeda identification grid 

As with the previous analysis method, the particle identification is per- 
formed using the AE-E method. In such a representation, the particles are 
grouped around quasi-hyperbolas characterized by their charge and mass. 
Then, for each detected particle, the closest fine has to be identified and the 
corresponding charge and mass have to be assigned to the particle. It can be 
done by using graphical contours as presented in the previous section or by 
using the KaliVeda identification tools. 

10 



KaliVeda is using identification grids which are graphs that contain a set 
of identification hnes as described by equation [2] (cf Fig. [sjand Fig. |9]). Each 
hue corresponds to an isotope and reproduces the shape of its quasi-hyperbola 
on the AE-E map. 
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Figure 8; Handmade identification grid built on a AE-E map for the most produced 
isotopes. 
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Figure 9: New identification grid obtained by fitting the Tassan-Got functional to the 
handmade identification grid. 

In a first step, to construct the grids, the user has to do a "handmade 
identification grid" for a few lines on the AE-E map. To do so, one has to 
give the coordinates of about ten points of the relevant lines and to specify 
charge and mass of the corresponding particle as shown on Fig. |8j Only few 
lines are needed here (those with the highest statistics), whereas the other 
one will be computed afterwards. 

In a setup with three stages telescopes, two identification grids are needed 
per telescope (one for the particles stopped in the thick Silicon detector 
and one for those stopped in the Csl). Fig. M represents the "handmade 
identification grid" for the particles stopped in the thick Silicon detector 
(superimposed to the AEthm-Ethick map). 

In a second step, the functional of equation |2] is fitted to the handmade 
AEthin-Ethick identification grids. The same method is used for the AEthick- 
Ecsi identification grids but the energy deposited in the Csl is here expressed 
as a function of the emitted light as described in equation |3} The fit result 
allows to obtain the analytical expression of the identification lines for all 
isotopes. New identification grids including the lines of isotopes with low 
statistics can now be created, as shown on Fig [9j 

By this identification method, two sets of accurate AE-E identification 
grids (AEthin-Ethick and AEthick-Ecsi) have been determined. 

12 



4-2. The KaliVeda method 

In this part, the method we have developed will be presented. This 
method is using the KaliVeda tools described above. 

4.2.1. Particle identification 

To identify the detected particles, the procedure works as follows: 

• For each experimental measurement, the associated identification grid 
is loaded and the algorithm computes the distance on the AE-E map 
between the coordinates of the detected particle ((AEthin;Ethick) or 
(AEthick;Ecsi)) and the different isotope's lines. 

• The charge and mass assigned to the particle are the one corresponding 
to the closest line of the identification grid. 

• A new variable named Areai is determined as a function of the distance 
to the isotope's line with A^ea/ G [A— 0.5;A+0.5]. The value of P^reai is 
close to the mass number A when the experimental point is close to the 
identification line. The minimal or maximal values of A^gai correspond 
to points located at the same distance of two lines with consecutive 
mass numbers. 

• This new variable is used to determine the Particle Identifier (PID) 
expressed as PID = Z + 0.1 x(Area/ - 2xZ). 

As examples, the PID for a particles may range from 1.95 to 2.05, the PID 
of ^He isotopes may ranges from 1.85 to 1.95 and the PID for ^He may range 



from 2.15 to 2.25. Fig. 10 illustrates an example of PID calculation for alpha 



particles. The result of the PID calculation for all the detected isotopes is 



shown on Fig. 11 The quality of the identification is very good and the 



different isotopes are clearly separated, even those with low statistics. 
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Figure 10: Example of PID calculation for ''He fragments. 




Figure 11: PID calculation for all detected isotopes. 



4-2.2. Energy calibration 

Once the particle identification has been achieved, the energy cahbration 
of the detectors can be performed. 

For energy cahbration of the Sihcon detectors, the analytical expres- 
sions (eq. [2| of the identification grids for the AEthm-Ethick map are used. 



14 



First, a theoretical identification grid is generated. The points of this grid 
are obtained with the energy loss calculation tool of KaliVeda, taking into 
account the Silicon detectors thicknesses. The energy loss functional is then 
fitted to the theoretical identification grid. It gives an analytical expression 
of the theoretical lines. Then, the analytical expression of the experimental 



lines are obtained as described in part 4.1.2 



Two sets of parameters have thus been obtained, one for the theoreti- 
cal grid expressed in MeV, and one for the experimental one expressed in 
channels. As Silicon detectors are delivering a linear response, their energy 
calibration consists in determining the afiine relation between channels and 
MeV. Four calibration parameters have to be determined, two for the thin 
Silicon detector and two for the thick Silicon detector. By applying the cor- 
rect transformations respectively to AE and E in the experimental energy 
loss functional, the two functionals (experimental and theoretical) will be 
equal. This procedure is done by using an iterative algorithm which adjusts 
the four calibration parameters to make the two functionals as close as possi- 



ble. Fig. 12 represents a AE-E calibrated map on which the theoretical grid 
is superimposed. Both are in perfect agreement. 
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Figure 12: Result of the Silicon detectors calibration. The theoretical grid is superimposed 
to the AE-E calibrated map (the right-hand side is a zoom of the red square). Only 
statistical errors arc represented. 



The Csl energy calibration can not be done with the same method because 
of the non linear response of the emitted light versus the deposited energy 
(cf eq. |3|. The thick Silicon detector energy calibration is used to determine 
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the deposited energy in the Csl as follows. Once a particle is identified in 
the AEthick-Ecsi rnap, its charge, mass and energy loss in the thick Silicon 
detector are known. Then, a KaliVeda tool computes the corresponding 
incident energy of the particle thanks to the known energy loss in the Silicon 
detector. The residual energy which will be deposited in the Csl can thus be 



deduced. The disadvantage of such a method, as the one used in section |3.2 
is that the error made on the thick Silicon detector calibration parameters is 
propagated to the energy deposited in the Csl. This error has been estimated, 
in this case, from around 20% for 95 MeV protons to 4% for 95 MeV/u ^^C. 

The main advantage of using the KaliVeda toolkit is that once the identi- 
fication grids are created, they allow the detector calibration and the particle 
identification at the same time. With such a method, the energy calibration 
is done not only on a few points but for all the points of the identification 
lines. Moreover, these procedures are automated which saves time in data 
analysis. The time required, for the user, to construct the necessary elements 
for the calibration and identification is about one hour per telescope. 

5. Comparison of the two analysis methods 

The obtained data have been analyzed with the two methods described 
above. In this part, some results of fragmentation cross section measurements 
for the carbon target will be presented and the two analysis methods will be 
compared. 

Firstly, the angular distributions, obtained for angles from 4° to 43°, will 
be compared. The ^^C fragmentation cross sections for a ^X fragment are 
obtained as follows : 



zX) = ., ,,o,, .,,A,, .r (4) 



da . _ Nax X Atarget 

dQ^"^ ' ~ Ni2c xQx pxthx A/'a 

where Nax is the number of ^X fragments detected, Atarget is the target 
mass, Ni2q is the number of incident carbon nuclei, Q is the solid angle of 
the detector, p and th are the target density and thickness respectively and 
A/'a is the Avogadro's number. 

Fig. 13, Fig. 14 and Fig. [15] display the angular distributions obtained with 



the two analysis methods for different emitted fragments. Statistical errors 
only are represented, systematic errors are still under study. As expected, the 
fragment emission is dominated by H and He isotopes with a predominance 
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of a at low angles (below ~ 10°) which is compatible with the alpha cluster 
structure of the ^^C The results also show an angular emission most forward 
peaked for heavier fragments. 
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Figure 13: Comparisons of the Hydrogen angular distributions obtained with the two 
analysis methods. Only statistical errors are represented. 
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Figure 14: Comparisons of the Helium angular distributions obtained with the two analysis 
methods. Only statistical errors are represented. 
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Figure 15: Comparisons of the angular distributions obtained with the two analysis meth- 
ods for the heavier fragments. The Z distributions are obtained by adding the cross 
sections of each isotope. 



The two analysis methods give consistent results. The discrepancies be- 
tween the two methods are of ± 5% for most of the isotopes and angles. 
However, the discrepancies can reach 15% at low angles for light isotopes 
or for isotopes with low statistics. These differences are due to the widths 
of the particle selection cuts which are not the same for the two methods 
(width of the contours for the first method and width of the identification 
lines for the second one). The differences are more important for protons be- 
cause of the increasing "background" on the AE-E maps at low energy. As 
this "background" is mainly due to protons escaping the Csl by scattering, 
a large identification line width has been used in the KaliVeda method in 
order to not lose proton statistics. 



Secondly, the energy distributions will be compared. Fig. [16] represents 
the energy distributions of protons and a fragments at different angles ob- 
tained with the two analysis methods. These distributions are dominated at 
low angles by a peak centered close to the beam energy (95 MeV/u). The 
energy and amplitude of this peak are decreasing with the emission angle. 
This observation confirms that most of the detected fragments are coming 
from the projectile fragmentation. 
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Figure 16: Comparisons of the energy distributions obtained with the two analysis methods 
for protons and a at different angles. 



The energy distributions are also in good agreement. The discrepancies 
between the two methods are due on one hand to the propagation of the 
energy cahbration of the thick Sihcon detector. This effect can be seen on 
Fig. 16(b) in which the differences between the two analysis method are 



increasing with the energy. It should be noticed that the discrepancies stay 
under 10% even for high energy protons. On the other hand, the discrepancies 
are also coming from the particle selection width as explained for the angular 
distribution. This effect is clearly visible for the higher energies of Fig. |16(a) 
in which the largest width of the KaliVeda methods leads to the selection of 
background (due to the presence of ^H and ^H). As the background of the 
experiment is low, this effect is statistically low even with the choice of a 
selection width larger in the KaliVeda method than in the classical one. The 
comparisons between the energy distributions obtained with the two analysis 
methods show the coherence and the robustness of the two identification and 
calibration methods. 



6. Conclusion 

Two different methods for the data analysis of impinging ^^C ions at 
95 MeV/u on thin carbon target have been presented in this paper. Mea- 
surements for twenty different angles led to production rates, angular distri- 
butions and energy spectra for fragmented particles, ranging from protons 
to carbon ions. Two different analysis methods have been used, the first one 
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is based on graphical cuts while the second one uses the KaliVeda frame- 
work. Both of them are in good agreement, within a discrepancy close to 5% 
for most angular distributions, the discrepancy is higher concerning energy 
distributions in some specific cases. Angular distributions are dominated by 
the emission of light fragments (Z < 3), more especially by a particles at 
forward angles, that is consistent with the a cluster structure of the ^^C ion. 
Energy distributions dominated at low angles by a peak close to the beam 
energy indicate that most of the emitted particles result from the projectile 
fragmentation during nuclear reactions. As mentioned above, other targets 
have been used in this experiment using the same beam at the same energy, 
which will be the subject of an upcoming paper. The analysis of these data 
will help to determine the double differential cross section for the different 
couples C-H, C-0, C-Al and C-Ti. 
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